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ABSTRACT 

We report the discovery of the new gravitational lens system CLASS B0631+519. 
Imaging with the VLA, MERLIN and the VLBA reveals a doubly-imaged flat- 
spectrum radio core, a doubly-imaged steep-spectrum radio lobe and possible 
quadruply-imaged emission from a second lobe. The maximum separation between 
the lensed images is 1.16 arcsec. High resolution mapping with the VLBA at 5 GHz 
resolves the most magnified image of the radio core into a number of sub-components 
spread across approximately 20 mas. No emission from the lensing galaxy or an odd 
image is detected down to 0.31 mJy (5a) at 8.4 GHz. Optical and near-infrared imag- 
ing with the ACS and NICMOS cameras on the HST show that there are two galaxies 
along the line of sight to the lensed source, as previously discovered by optical spec- 
troscopy. We find that the foreground galaxy at z=0.0896 is a small irregular, and 
that the other, at z=0.6196 is a massive elliptical which appears to contribute the 
majority of the lensing effect. The host galaxy of the lensed source is detected in 
the HST near-infrared imaging as a set of arcs, which form a nearly complete Einstein 
ring. Mass modelling using non-parametric techniques can reproduce the near-infrared 
observations and indicates that the small irregular galaxy has a (localised) effect on 
the flux density distribution in the Einstein ring at the 5-10% level. 

Key words: gravitational lensing; cosmology. 



1 INTRODUCTION 

Strong gravitational lensing occurs when a background 
source is multiply imaged by the gravitational field of an 
intervening massive galaxy or cluster of galaxies. Observa- 
tions of gravitational lenses have been used in many astro- 
physical applications, such as in modelling the mass dis- 
tributions of di stant galaxies and as constraints on cosmo - 
logical models llKochanek. Schneider fc Wambseans'sl|2004) 
CLASS (Cosmic L ens All-Sky Survey; iMvers et alJ l20oi 
iBrowne et all 120031) and JVAS (Jodr e ll Bank - VLA As- 
trometric Survey; |Patnaik et al]ll992t IBrowne et al J Tl*998t 
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I Wilkinson et all Il998t) are radio surveys of the northern 
sky designed to find gravitationally lensed compact radio 
sources. In JVAS and CLASS 11685 sources in a statisti- 
cally well-defined sample have been examined and 22 lens 
systems found. In this paper, we announce the discovery 
of the last of these systems, CLASS B0631+519. This lens 
system exhibits complex radio structure over scales from 
3.6 mas to 1.16 arcsec and has a nearly co mplete inf rared 
Einstein ring, similar to tha t in B 1938+666 jPatnaik et alJ 
11991 iKing et alJll997l Il998t) . It is a system with one of the 
richest lensed image structures known and is thus an ideal 
system with which to probe mass properties of the lensing 
galaxy. B0631+519 is a member of the CLASS statistically 
complete sample. 
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In Section [5] we discuss radio observations of B0631+519 
made using the Very Large Array (VLA), Very Long Base- 
line Array (VLBA) and the Multi- Element Radio-Linked In- 
terferometer Network (MERLIN; |Thomassor]|l986l) . Section 
[3] presents optical and near- infrared images taken with the 
William Herschel Telescope (WHT) and the Hubble Space 
Telescope (HST). In Section [| we consider the data in the 
context of the lens hypothesis and test some simple lens 
models, and in Section we conclude by discussing the fu- 
ture priorities for work on this lens. 



2 RADIO OBSERVATIONS 



As d escribed in IMvers et alJ (2003) and IBrowne et all 
(2003), potential lens candidates from CLASS are tested by 
observing them with a range of angular resolutions and at 
several different frequencies. If the supposed multiple images 
in a lens candidate are found to have very different surface 
brightnesses at any angular resolution the lens hypothesis 
can usually be rejected immediately, since lensing preserves 
surface brightness. In CLASS, candidates discovered with 
the VLA are observed with higher resolution using MER- 
LIN, and, if not eliminated at this stage, are re-observed 
with even higher resolution using the VLBA. Accordingly, 
radio observations of B0631+519 have been made over fre- 
quencies ranging from 1.7 GHz to 15 GHz, with angular 
resolutions from 3.6 mas to 236 mas. The observations are 
summarised in Table Q 



2.1 VLA 

B0631+519 was observed at 8.4 GHz on 1994 March 05 and 
1999 August 03. The 1999 observation was a follow-up to the 
1994 snapshot observation and had almost 5 times its sensi- 
tivity (Table 0. The snapshot was made using a bandwith 
of 50 MHz, and 3C48 was used to set the flux density scale; 
the follow-up observation used a bandwidth of 100 MHz, 
and the flux density scale was set using 3C147. 

All VLA observations of B0631+519 were made with the 
VLA in 'A' configuration, and were calibrated using AIPS 
and map ped with the C altech difference mapping program 
DIFMAP llShepherdl ll997 > ) . Model component positions and 
flux densities from the VLA observations are displayed in 
Tabled 

The deconvolved map from the 1999 follow-up observation is 
shown in Figure We label the detected radio components 
as shown in Figure Components Al and B are separated 
by 1.15 arcsec at a position angle (from brightest to dimmest 
source) of +135°, and the A1:B flux density ratio is 8:1. 
Components Al, A2 and B are unresolved. Component X is 
a knot of elevated surface brightness embedded in the jet- 
or arc-like feature that branches to the south of A2. There 
is some extended emission visible on the opposite side of 
Al from X, which we label Y, although the detection is 
marginal (about 2.3<r). The total flux density in the model 
is 42±2 mjy (where we have assumed 5% errors on absolute 
flux density measurements from the VLA at 8.4 GHz). 



AIPS: NRAO's Astronomical Image Processing System 
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Figure 1. Radio map of the VLA 8.4 GHz data taken on 1999 
August 03. The RMS noise in the residual map is 86 fijy beam -1 . 
The data are naturally weighted. The synthesized beam is sized 
250x223 mas with the major axis oriented along a position an- 
gle (P.A.) of —17.5°. The lowest contours are plotted at -3 and 
3 times the RMS noise level; successive contours represent an 
multiplicative increase in surface brightness by a factor of 1.5, 
up to the highest contour level of 22.8 mjy beam -1 . Compo- 
nent Al is located at RA: 06 h 35 m 12.31402 s , Dec: +51° 57' 
01.8026" (J2000.0) as measured by the VLBA. 



In the snapshot Al, A2 and B were detected, while X and 
Y were not. The snapshot shows an ALB flux density ratio 
of about 7:1, and a total flux density from the candidate of 
46±2 mjy. This is 10% higher than the equivalent quantity 
in the follow-up map. The discrepancy might be accounted 
for by uncertainty in the absolute amplitude calibration be- 
tween the 1994 and the 1999 observations, but it could also 
indicate that the source is variable. The mean size of the 
synthesized beam varies by less than 3 mas between the two 
(uniformly weighted) maps, or by 2% in area, so it is unlikely 
that 10% of the emission has been resolved out in one map 
with respect to the other. The ALB flux density ratio also 
changes by 16% between the two observations, which is hard 
to reconcile with a global mis-calibration of the amplitude 
scale. The signal-to-noise ratios on Al and B were at least 
10 in all VLA observations. 

A further observation of B0631+519 was made with the VLA 
at 15 GHz on 1999 July 02 for 354 seconds, with a bandwidth 
of 100 MHz. The amplitude scale was set with observations 
of 3C286. The deconvolved, naturally weighted map is shown 
in Figure [5] The total flux density detected at 15 GHz is 
34±3 mjy, and the ALB flux density ratio is 7.5:1. Both Al 
and B are unresolved by the VLA at this frequency. A2 is 
not detected, suggesting that it has a steeper spectral index 
than Al and B. 
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Date observed 


Instrument 


Frequency [GHz] 


Integration time 


Map residual noise [/^Jy beam 1 ] 


Mean beam size [mas] 


1994 March 05 


VLA 


8.4 


16 sec 


410 


234 


1999 July 02 


VLA 


15 


354 sec 


360 


112 


1999 August 03 


VLA 


8.4 


510 sec 


86 


236 


1999 December 10 


VLBA 


1.7 


1.4 hr 


130 


8 


2001 March 21 


VLBA 


5 


1.7 hr 


82 


3.6 


2001 July 07 


MERLIN 


5 


19 hr 


120 


52 


2002 December 01 


VLBA 


1.7 


7.8 hr 


64 


9 


2003 March 01 


MERLIN 


1.7 


20 hr 


62 


173 



Table 1. Radio observations made of B0631+519 to date. 



Component 


Frequency [GHz] 


Position offsets [mas] 
AR.A. ADec. 


Flux density [mjy] 


Al 


8.4 


+0±1 


+0±1 


34.3±1.7 




15 


+0±2 


+0±2 


29.6±3 


A2 


8.4 


-297±10 


-300±10 


2.2±0.1 




15 






<1.8 


B 


8.4 


+819±5 


-819±5 


4.2±0.2 




15 


+813±19 


-805±19 


3.95±0.4 


X 


8.4 


-114±33 


-754±33 


0.98±0.1 




15 






<1.8 


Y 


8.4 


+577±70 


+81±70 


0.2±0.1 




15 






<1.8 



Table 2. Model components from radio observations of B0631+519 with the VLA. All position offsets are relative to Al. The 8.4 GHz 
flux densities and positions are derived from the 1999 follow-up data, not the 1994 discovery data. Non-detections are quoted at the 5<r 
level throughout this paper. Errors on positions of Gaussian model components are estimated by dividing the beam size by the peak 
signal-to-noise ratio (SNR). Errors on flux densities are assumed to be subject to a 5% random error in absolute flux density calibration 
at 8.4 GHz and a 10% error at 15 GHz. 



2.2 MERLIN 

B0631+519 was observed at 5 GHz using MERLIN on 2001 
July 07, over a period of 19 hours. Telescope gains were 
calibrated using the source OQ208, and the observations 
were phase referenced to the source JVAS J0642+5247. All 
MERLIN data were reduced by initial amplitude calibra- 
tion and flagging with the MERLIN D-programs, followed 
by phase and amplitude calibration and further flagging in 
MPS. difmap was used to map and self-calibrate the data. 
3C286 was used to set the flux density scale for all MERLIN 
observations. Model component positions and flux densities 
from the MERLIN observations are shown in Table [3] 

The deconvolved 5 GHz map is shown in Figure |3] Com- 
ponent B has been resolved into two components (Bl and 
B2). Al and Bl appear compact, but A2 and B2 are slightly 
resolved. X and Y have been resolved out and there is no 
sign of the extended emission revealed by the 8.4 GHz VLA 
observations. 

Additional MERLIN observations were made at 1.7 GHz 
on 2003 March 01. Phase referencing was carried out us- 
ing the source J0636+5009. B0552+398 was used to cal- 
ibrate the telescope gains. The newly upgraded 76 metre 
Lovell Telescope was included in the interferometer net- 
work, allowing an RMS noise level in the deconvolved map of 
62 fj, Jy beam -1 to be reached in 20 hours. The deconvolved 
1.7 GHz map is shown in Figure 0] The map shows emission 
both to the west and east of image Al. We identify these 
regions of emission as low-frequency counterparts to compo- 



Al 



B 
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Figure 2. 15 GHz VLA data taken on 1999 July 02. The synthe- 
sized beam is 152x132 mas oriented with the major axis at a P. A. 
of +26°, and the RMS noise in the residuals is 360 fj,Jj beam -1 . 
The data are naturally weighted. Contours are plotted as for Fig- 
ure 1; the highest contour value is 27.8 mjy beam -1 . 
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Figure 3. Radio map of the 5 GHz MERLIN observations taken 
on 2001 July 07. The data were naturally weighted. The beam 
size is 58x46 mas with the beam's major axis being oriented at a 
P.A. of -34°. The RMS noise in the residuals is 120 /ijy beam -1 . 
Contours are plotted at levels of (-3, 3, 4.5, 6.75, 10.1, 15.2, 22.8, 
34.2, 51.3, 76.9, 115, 173, 259) times the RMS noise. 

nents X and Y seen in the VLA data. Bl also develops an 
extension to the west, which we label Z. 

2.3 VLBA 

The VLBA was used to observe B0631+519 at 5 GHz on 
2001 March 21 over a period of 1.7 hours. Data were taken 
using four contiguous 8 MHz bands and a single hand of cir- 
cular polarisation. The observations were phase-referenced 
using the calibrator source J0631+5311 and employed 2 bit 
sampling during digitisation; correlation was performed us- 
ing the VLBA correlator in Socorro. During correlation each 
band was subdivided into 16 0.5 MHz channels and the data 
were averaged into 2 s integrations. 

All VLBA observations were reduced using AIPS to flag 
and calibrate the data, following standard recipes for the 
VLBA. The aips task scim g was then used to map the data 
through cycles of CLEAN dHogbomlll974h followed by self- 
calibration. Flux densities and positions of VLBA compo- 
nents are shown in Table 01 

The deconvolved 5 GHz maps are shown in Figure [5] Uni- 
form weighting was employed. The resolution achieved was 
3.6 mas, sufficient to split Al into four sub-components (la- 
belled Ala to Aid) spread over an angle of about 20 mas. 
Bl remained unresolved, and no emission from the A2, B2, 
X, Y or Z components was detected down to the 5cr limit 
of 390 fiJy beam -1 , indicating that these components have 
probably been resolved out. 

We note that the sum of the flux densities of components 
Ala to Aid (46±2 mjy) is consistent with the MERLIN 
5 GHz flux of component Al to within the measurement er- 
rors, which implies both that Al did not vary by more than 



Al 
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Figure 4. Radio map of the 1.7 GHz MERLIN observations taken 
on 2003 March 01. The data were naturally weighted. The syn- 
thesized beam is sized 186x161 mas at P.A. —31.7°. Contours are 
plotted at factors of (-3, 3, 4.5, 6.75, 10.1, 15.2, 22.8, 34.2, 51.3, 
76.9, 115, 173, 259, 389, 584, 876) times the RMS noise. 

5% over the 4 months separating the observations and also 
that no radio structure seen in Al by MERLIN was resolved 
out by the VLBA, despite the factor of ~14 increase in res- 
olution. However, another possible explanation is that an 
increase in source flux density (due to intrinsic variability) 
cancelled out any loss of emission due to resolution effects. 
The same conclusions and caveats apply to the MERLIN 
and VLBA observations of Bl at 5 GHz. 

Comparison of the total Al flux density to the flux density 
of Bl suggests an A1:B1 magnification ratio of ~9:1, and 
hence that sub-components in Bl corresponding to those in 
Al would span about 2 mas. Resolving the expected sub- 
components in Bl is thus well within the capabilities of 
VLBI (Very Long Baseline Interferometry) . 

Two epochs of phase-referenced VLBA data were obtained 
at 1.7 GHz, on 1999 December 10 and 2002 December 01. 
The phase reference source used was J0631+5311. The 1999 
data was taken in two polarisations via 4 contiguous bands, 
2 bands to each polarisation. 2-bit sampling was used, and 
correlation at Socorro again subdivided each band into 16 
0.5 MHz wide channels and averaged the data into 2 s inte- 
grations. The phase centre was taken halfway between im- 
ages Al and Bl to minimise bandwidth smearing. 

The CLEANed, naturally weighted maps of Al, A2 and 
B1/B2 from the 2002 epoch are shown in Figure|S| Al ap- 
pears resolved but (unsurprisingly) smoothed compared to 
the 5 GHz VLBA map, while Bl remains unresolved. A2 
is resolved into an arc-like feature of low surface brightness 
that extends over 80 mas. B2 is partially resolved into an ex- 
tended source. No emission was detected from components 
X, Y or Z. 

The maps from the two 1.7 GHz epochs look very similar. 
The flux densities of Al, A2 and B2 agree between the two 
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Component Frequency [GHz] Position offsets [mas] Flux density [mjy] 

AR.A. ADec. 



Al 


1.7 


+0±1 


+0±1 


66.4±6.6 




5 


+0±1 


+0±1 


46.9±4.7 


A2 


1.7 


— 293±1 


— 313±1 


19±1.9 




5 


— 293±3 


— 308±3 


4.2±0.4 


B 


1.7 


+809±1 


-800±1 


11.8±1.2 


Bl 


5 


+828±1 


-816±1 


5.4±0.5 


132 


5 


+801±8 


-722±8 


1.3±0.1 


X 


1.7 


-40 


-880 


11±1 




5 






<0.6 


Y 


1.7 


+595±6 


+67±6 


2.7±0.3 




5 






<0.6 


Z 


1.7 

5 


+691±5 


-1016±5 


5.0±0.5 
<0.6 



Table 3. Radio components detected in B0631+519 with MERLIN. All position offsets are relative to Al. X is a large arc that is not well 
described by a single gaussian component; the position quoted is an estimate for the centre of the feature. The integrated flux density of 
X was derived by summing within the 3<r contour. We have assumed a 10% absolute calibration error on flux densities from MERLIN. 
Position errors are estimated from the beam size divided by the peak SNR. 



epochs to within the measurement errors. The flux density 
of Bl decreases by about 9% between 1999 and 2002; this 
decrease is not quite significant at the 2a level given the 
measurement errors of about 5% for Bl's flux density. The 
measurements of the Al/Bl image separation from the two 
epochs are consistent to within 0.5 mas, while the A1-B2 
separation measurements agree to within 2 mas. 

To identify which of the four VLB A sub-components of Al 
correspond to the flat-spectrum radio core in the lensed 
source, we registered the 1.7 GHz map (from the 2002 epoch) 
with the 5 GHz map by using the position of component Bl 
as a reference (since Bl is largely unresolved in both maps). 
To overlay Bl in this way required the 1.7 GHz image to 
be translated by 1.7 mas south and 8.9 mas east because 
slightly different positions were used for the phase reference 
source J0631+5311 between the 5 GHz and 1.7 GHz obser- 
vations. The translated 1.7 GHz map is shown overlaid on 
the 5 GHz map in Figure [7| Sub-component Ala appears 
to be the flat spectrum core, while the peak in the 1.7 GHz 
map corresponds to component Ale. Therefore at 5 GHz we 
apparently see both approaching and receding jets in the 
source AGN. 



2.4 Radio spectra from continuum data 

We plot the radio spectra for the model components as a 
function of frequency in Figure |H| If the lens hypothesis is 
correct then all the lensed images of a source should have 
similar spectra, provided that the lens galaxy ISM does not 
scatter or absorb th e source light, and that t he source flux 
density is constant jKochanek fc Dalall Eool . Al and Bl 
have similar spectra, as do A2 and B2. Possibly X, Y and Z 
also have similar spectra, although at 8.4 GHz component Y 
is a marginal detection and component Z is a non-detection 
(we have plotted a 5a upper limit). We suggest that there 
are three separate source regions revealed by the radio obser- 
vations, two of which are doubly-imaged and one of which is 
quadruply imaged (generating images Y, Z and the arc X). 
We will discuss the geometry of the source with respect to 
the lens caustics further in Section 2] 



3 OPTICAL OBSERVATIONS 

Typically lens galaxies are radio-quiet ellipticals, so op- 
tical/infrared imaging is necessary to establish the lens 
galaxy position when the lenses are identified through ra- 
dio searches. The position of the lens galaxy with respect 
to the lensed images can be used to constrain mass mod- 
els, and in some systems unique information on the lensed 
source is available through optical imaging, for instance 
where an optical or infr ared Einstein ring is vis ible as seen 
in CLASS B1938+666 jKing et aljll997l. Il99ctl . We have 
observed B0631+519 at optical/near-infrared wavelengths 
with the WHT through an R-band filter and with the HST 
through F555W, F814W and F160W filters. Details of these 
observations are listed in Tabled 

3.1 WHT observations 

B0631+519 was observed with the 4.2 metre William Her- 
schel Telescope (WHT) on 2002 February 03, through an ri- 
band filter. The detector used was a 1024x1024 pixel TEK 
CCD, having a total field of view of 111 x 111 arcsec. The 
data were calibrated using the Starhnk package CCDPACK. 
B0631+519 was partially resolved, but there was insuffi- 
cient resolution to separate any lensed images from the lens 
galaxy. 

Aperture photometry was undertaken on the WHT image 
using ob servations of L andolt standard stars 93 407 and 
100 280 jLandoltlliggl to correct for atmospheric extinc- 
tion and to calibrate the brightness scale in magnitudes. An 
elliptical aperture with major diameter 4.3 arcsec and an 
axis ratio of 0.75 was used to obtain an R-band magnitude 
for B0631+519 of 21.1 ±0.1, corrected for Galactic extinctio n 
using the dust maps of lSchleeel. Finkbeiner fc Davisl |l998). 

3.2 HST optical and near-infrared observations 

B0631+519 was observed with the Hubble Space Telescope 
(HST) on 2003 August 19 and 2004 February 15 as part of 
HST proposal 9744, "HST Imaging of Gravitational Lenses" 
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Component Frequency [GHz] Position offsets [mas] Flux density [mjy] 

AR.A. ADec. 



Al 


1.7 


+0±0.1 


+0±0.1 


51.4±3 


Ala 


5 


+0±0.1 


+0±0.1 


20.5±1 


Alb 


5 


+5.7±0.1 


+3.5±0.1 


8.8±0.4 


Ale 


5 


+11.1±0.1 


+6.6±0.1 


7.3±0.4 


Aid 


5 


-7.2±0.1 


-5.0±0.1 


9.2±0.5 


A2 


1.7 


-300 


-306 


6.6±0.7 




5 






<0.39 


Bl 


1.7 


+812.5±0.1 


-827.5±0.1 


5.6±0.3 




5 


+822.5±0.1 


-821.5±0.1 


5.1±0.3 


B2 


1.7 


+790.7±0.4 


-732.9±0.5 


2.2±0.2 




5 






<0.39 



Table 4. Model components for the radio observations of B0631+519 with the VLBA. Stated errors on flux densities assume a 5% 
random error in the absolute flux density scale calibration at the VLBA. The 1.7 GHz position offsets are given relative to the surface 
brightness peak in Al. For A2 and B2 at 1.7 GHz, we have doubled the error to reflect the uncertainty in the flux densities of these 
resolved components; their flux densities were obtained by summing within the 3<r contour. Positions for Al, Bl and B2 were determined 
from the best-fit gaussians, restricted to the compact north-eastern feature in the case of Al at 1.7 GHz. The position given for A2 is 
that of the brightest pixel. For the 5 GHz data, position offsets arc taken relative to Ala at RA: 06 h 35 m 12.31390 s , Dec: +51° 57' 
01.7949" (J2000.0). 



Instrument 


Detector 


Date observed 


Filter 


Integration time [sec] 


PSF FWHM [arcsec] 


Plate scale [arcsec pixel 1 ] 


HST 


ACS/WFC 


2003 August 19 


F555W 


2236 


0.05 


0.05 


HST 


ACS/WFC 


2003 August 19 


F814W 


2446 


0.08 


0.05 


HST 


NICMOS/NIC2 


2004 February 15 


F160W 


2560 


0.13 


0.075 


WHT 


Aux. port CCD 


2002 February 03 


R 


1800 


0.7-0.9 


0.108 



Table 5. Optical and infrared observations made of B0631+519 to date. 



(PI: C.S. Kochanek). Images were taken using the Advanced 
Camera for Surveys (ACS) Wide Field Channel (WFC) 
through the F555W and F814W filters. B0631+519 was 
also observed with the NICMOS camera NIC2 through the 
F160W filter, in the MULTIACCUM observing mode. The 
wideband HST filters F555W, F814W and F160W approx- 
imately correspond to Johnson V, I and H respectively. 
Details of the optical observations are shown in Table 
The ACS observations used a dither pattern consisting of 
two pointings separated by 0.73 arcsec. Two exposures were 
taken at both pointings to simplify cosmic ray rejection. The 
NICMOS data used a four-point dither pattern that placed 
B0631+519 at the centre of each quadrant of NIC2 in turn. 

The ACS and NICMOS data were debiased, dark- 
subtracted and flat-fielded using the OTFR (On-The-Fly- 
Recalibration) facility of the MAST HST archive. W e then 
used the multidrizzle script (iKoekemoer et all 12002') from 
within PyRAF to correct the ACS data for distortion, to 
eliminate cosmic rays and to combine dithered pointings 
into a single output frame, the shifts between frames be- 
ing obtained from application of cross-correlation through 
the crossdriz task in IRAF. The NICMOS data were com- 
bined using pydrizzle to interpolate the data onto an output 
image with a uniform pixel scale of 0.075 arcsec, using the 
standard distortion solution for NIC2. Cosmic rays in the 
NICMOS data were identified and blanked by the OTFR 
pipeline itself. The IRAF task pedsub was used to remove 
some residual dark signal that was evident in the output 
from the pipeline. Both the ACS and NICMOS images were 



aligned to the local celestial coordinate axes as part of the 
drizzling process. 

The HST images (Figure |5J reveal two galaxies at the lens 
system position, surrounded by arcs of emission from the 
lensed source at longer wavelengths. We label these galaxies 
LI and L2. LI appears to be an elliptical galaxy, and is de- 
tected in all three bands. In the F814W and F160W images 
some emission from the source appears in the form of arcs 
of extended emission encircling LI, consistent with the lens 
hypothesis if LI is the primary lensing galaxy. L2 appears 
to be an irregular galaxy and is seen in the F555W and 
F814W data, although in the latter it is confused with one 
of the lensed arcs. L2 is not cleanly detected in the NICMOS 
image, and may not be visible in the near-infrared. 

The ACS images have a large enough field of view to show 
the environment of B0631+519. The nearest extended object 
is located 4.4 arcsec away at a P. A. of -72°, and has a total 
Imug ~ 22.3 in a 4 arcsec diameter aperture. The object's 
light distribution is not smooth and is difficult to classify. 

We carried o ut surface-br ightness fitting on LI using a Sersic 
light profile JSersidll968j)_ convolve d with appropriate PSFs 
produced by Tiny Tim (|Krjsg|l993j) . The profile parameters 
were varied to minimise the sum of squared residuals be- 
tween the data and the model. When fitting Li's surface- 
brightness profile we masked out L2 and the lensed images 
to avoid disturbing the fit. L2 and the lensed arcs do not 
have easily parametrizable surface light distributions, so we 
did not attempt to separate them in the F814W or F160W 
images, and report only the sum of light remaining after 
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Figure 9. HST images of B0631+519. The top row of images show the lens before subtraction of Ll, while the bottom row shows the 
residuals after subtraction of the best- fit Sersic profile. North is up and east is left. The images are 3x3 arcsec in size, and are centred 
on Ll. 



Parameter Filter 

F555W (V) F814W (I) F160W (H) 



Offset in R.A. [mas] +360±80 +380±80 +360±80 

Offset in Dec. [mas] -500±80 -530±80 -500±80 

Total magnitude 22.7 20.0 17.7 

Effective radius [mas] 540 560 550 

Ellipticity 0.09 0.17 0.12 

Position angle -48 -53 -57 

Power-law index 0.29 0.28 0.30 



Table 6. Parameters of the best-fit Sersic profiles fitted to Ll in 
the three HST images. A power-law index of 0.25 corresponds to a 
de Vaucouleurs surface-brightness profile. The total magnitudes 
are specified in the Johnson photometric system, corrected for 
galactic extinction. Position offsets are given relative to image 
Al. 



Source 


V 


Brightness 
Ft I 


H 


Total 


22.1±0.1 


21.1±0.1 19.9±0.1 


17.3±0.1 


Ll 


22.8±0.1 


20.0±0.1 


17.8±0.1 


L2 


23.0±0.1 




>23 


L2+images 




22.3±0.1 




Images 


>25 




18.3±0.1 



Table 7. Results of photometry on B0631+519. The figures re- 
ported in this table were measured within a circular aperture of 
4 arcsec diameter. For the total luminosity of the subtracted Sersic 
models see Table© The WHT R-band magnitude was calculated 
using a Kron elliptical aperture 4.3 arcsec in (major axis) diame- 
ter. All figures have been corrected for galactic and, in the case of 
the WHT, atmospheric extinction. Non-detections are reported 
as lower limits on magnitude at the 5<r level. 



subtraction of Ll. In the F555W image the arcs were not 
detected, so we used GAIA to fit an elliptical aperture to 
L2 (after Ll had been subtracted). The centroid of L2 was 
found to be —0.33 arcsec east and +0.74 arcsec north of Ll. 
The axis ratio of the best-fit aperture was 0.56, with the ma- 
jor axis having a position angle —1.1° east of north. Table HJ 
shows the properties of the Sersic profiles subtracted from 
Ll. 

The ACS and NICMOS images following subtraction of Ll 
are shown in Figure^ The Sersic profile is a reasonable fit to 
the HST data, leaving maximum residuals of 3% and 4% of 
Li's peak brightness in the case of the F814W and F160W 
images respectively. The residuals are undetectable within 
the noise in the case of the F555W image, which has the 
lowest signal-to-noise ratio of all the HST observations. 

We carried out aperture photometry on B0631+519, using 
a circular aperture 4 arcsec in diameter to measure the flux 
of the lens system as a whole, of the Sersic model for Ll and 
of whatever light remained after subtraction of the Sersic 



model from the data. The measurements were converted into 
the Johnson magnitude system using published zero points 
for HST magn itude systems. Galactic ext i nction was esti- 
mated from ISchleeel. Finkbeiner fc Davis! <ll998l) . The re- 
sults are shown in Table Q The differences in Ll magni- 
tudes between Tables [5] and [7| are due to the fact that the 
values given in the latter are sums of the flux within finite 
apertures, while those in the former are equivalent to sums 
within infinite apertures. 

The lens galaxy position with respect to the lensed images, 
when known, provides two constraints on simple paramet- 
ric lens models under the assumption that the centre of 
(smoothly distributed) light and the centre of mass for a 
galaxy are coincident. Since the lens galaxy in B0631+519 
is apparently radio-quiet, the radio and HST data must be 
overlaid as accurately as possible to establish the lens posi- 
tion. The lensed images at near-infrared and optical wave- 
lengths exhibit no obvious point-like components that could 
be unambigiously aligned with the radio structure. We must 
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Figure 5. Radio map of the 5 GHz VLBA observations taken on 
2001 March 21, showing image A (top) and image B (bottom). 
The data were uniformly weighted. The beam-size is 3.8x3.4 mas 
oriented at P.A. +144° and the RMS noise is 82 /ijy beam -1 . 
Contours are plotted at levels of (-3, 3, 4.5, 6.75, 10.1, 15.2, 22.8, 
34.2, 51.3) times the RMS noise. 



therefore rely on absolute astrometric calibration based on 
stars visible in the field to align the optical and near-infrared 
images with radio data. 

We used the 2MASS 2 catalogue with Starhnk's gaia to 
establish an astrometric calibration for the ACS images 
(F555W and F814W filters). 2MASS positions have no sig- 
nifican t systematic offsets with respect to the radio-based 
ICRF Jda Silva Neto et al]f2005» . The astrometric solutions 
were based on twelve star positions and had reported RMS 
errors of ±0.08 arcsec in R.A. and declination, although ex- 
amination of objects around the field indicated positional 
agreement between the two calibrated images to within 
0.05 arcsec. In the case of NICMOS the relatively small field 
of view (19.5x19.5 arcsec) meant that only one star near to 
B0631+519 was imaged. The centre of this star was used as 
a reference to fix the absolute astrometry for the image. We 
estimate the astrometric error in the resulting solution for 



"i 1 1 r 








Bl 



Figure 6. Radio maps of the 1.7 GHz VLBA observations taken 
on 2002 December 01. The data were naturally weighted. The 
synthesized beam is an elliptical gaussian sized 11.2x7.5 mas ori- 
ented with its major axis at a P.A. of +177°. The RMS noise is 
64 fijy beam -1 . Contours are plotted at levels of (-3, 3, 4.5, 6.75, 
10.1, 15.2, 22.8, 34.2, 51.3, 76.9, 115, 173) times the RMS noise. 



Sec: http:/ /pegasus. astro. uuiass.edu 
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Figure 7. The 1.7 GHz VLBA observations from 2002 December 
01 plotted with the 5 GHz data from the same instrument. The 
5 GHz data are shown in greyscale, while the 1.7 GHz data are 
shown contoured (with the same contour levels as in Figure 151. 
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Figure 8. Plot of the radio spectra of the various radio compo- 
nents identified in B0631+519. The total flux density is based o n 
measurements at 327 M Hz from WENSS fRengelink et al.ll997^ . 
at 1.4 GHz from NVSS ICondon et a.lj|l99gl) and at 5 GHz from 
the GB6 survey jGregorv^^alT^99^K 



the drizzled NICMOS image at 0.08 arcsec, since the star's 
position on the sky was derived from the astrometrically 
calibrated ACS F555W image. The rotation of the image 
was checked by calculating the position angle of the vector 
from the field star to Li's position. This angle agreed with 
the same quantity measured from the ACS images to within 
0.05°. Figure ITTjl shows the 1.7 GHz MERLIN map overlaid 
on the ACS and NICMOS images. 

Fitting Sersic profiles to LI gave us best-fit pixel positions 
for the galaxy centre in all three HST filters, which we con- 
verted into celestial coordinates using the absolute astro- 
metric solutions from the HST images. The sky position of 
image Al (from the 1.7 GHz VLBA observations) was then 
subtracted to leave Li's offset on the sky from Al. These off- 



sets are reported in Table|S| The agreement between the off- 
sets from the F555W and F160W images (which effectively 
share the same astrometric solution) and the disagreement 
between F555W and F814W suggests that the error in the 
galaxy position is dominated by the random error in the 
astrometric solutions for F555W and F814W; error contri- 
butions from the galaxy subtraction procedure and from the 
distortion solutions for each camera are negligible by com- 
parison. The mean offset from Al to LI (based on the ACS 
data only) is +0.37T0.08 arcsec in R. A. and -0.52T0.08 arc- 
sec in declination. 



3.3 Redshift-dependent properties of B0631+519 

iMcKean et alJ i2004h have obtained an optical spectrum of 
B0631+519 with the W. M. Keck telescope. They report the 
discovery of two distinct galaxies at two different redshifts. 
The first is an emission line galaxy at z=0.0896±0.0001; 
and the second, which is the dominant source of emission 
in the spectrum, is an elliptical galaxy at z=0.6196±0.0004. 
These spectral classifications are consistent with the mor- 
phologies, colours and flux densities of the two foreground 
galaxies detected in the HST imaging. Therefore, we asso- 
ciate the 2=0.0896 emission line galaxy wi t h L2, and the 
2=0.6196 elliptical with LI. IMcKean etafl J2004T) refer to 
these galaxies in redshift order as Gl and G2 respectively; 
in this paper we have swapped the name order because the 
HST imaging demonstrates that LI (G2) is the primary lens- 
ing galaxy, something that was not obvious from the spec- 
troscopy alone. 

IMcKean et"afl J2004I) also discuss the tentative detection of 
a single emission line from the lensed source at 5040 A. As- 
suming that the line is either due to Lya or Mg II emission, 
they suggest 3.14 or 0.80 for the source redshift, respectively. 
However, further spectroscopy will need to be carried out 
to test these possibilities. We do not have enough (unam- 
biguous) photometric information on the lensed images to 
attempt to choose between these possibilites by using pho- 
tometric redshift techniques, but, in the absence of further 
spectroscopy, we can attempt to constrain the source red- 
shift by determining the rest-frame luminosity L of LI and 
converting it into a velocity disp ersion a using the funda- 
mental plane of ellip tical galaxies llDiorgovski fc DavisjigS^ : 
Dressler^J^iL^)87 (.This approach has been used before by 
Tonrv fc Kochanek i2000l) to estimate the source redshift in 



CLASS B1938+666. 

We converted the apparent magnitudes of LI and L2 as mea- 
sured by the HST (after translation into Johnson filters V, 
I and H, and following application of extinction corrections) 
into absolute magnitudes measured in Johnson B using 



M L 



mg - 5 log 



D 



Wpc 



- Kbq(z), 



(1) 



where Mb is the predicted absolute magnitude of the object 
in B, rriQ is the apparent magnitude of the object in some 
filter Q, Dl{z) is the luminosity distance 3 to the object at 
redshift z and Kbq(z) is the K-correction from the object's 



3 All cosmological distances in this paper were calculated assum- 
ing a flat universe with fl m = 0.27 and Q\ = 0.73. 
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Figure 10. The MERLIN 1.7 GHz radio data (contoured) over- 
laid on the NICMOS and ACS images (greyscale). The ACS im- 
ages cover 2.45x2.45 arcsec, while the NICMOS image covers 
2.475x2.475 arcsec. All images are centred on LI. North is up 
and east is left. 



redshift an d filter Q to the rest-f rame in B. We adopt the 
redshifts of iMcKean et al.l J2004I) for LI and L2. The K- 
corrections were calculated using the SYNPHOT package in 
IRA F, together wi t h tem plate spectra for LI and L2 from 
iBruzual fc Charlod <199Sl) and lKinnev et alJ (Il996l) . 

The K-correction for LI was obtained from the model ellipti- 
cal galaxy spectra of lBruzual fc Charlod jl99ot) . who provide 
multiple templates with various star formation histories and 
initial mass functions (IMFs) . We selected the template that 
produced the lowest scatter in the predicted rest-frame B 
magnitudes for LI amongst the three filters (F555W, F814W 
and F160W). The optimum template gave an RMS scatter 
amongst the three filters of ~ 0.02 mag, and was generated 
assuming a Salpeter IMF with mass range from 0.1 to 125 
Mq, in which star formation took place >5 Gyr ago during 
a burst lasting 1 Gyr. The templates were used solely to 
obtain a spectral K-correction; no correction for luminosity 
evolution was applied. 

The K-correction for L2 was e stimated using the st arburst 
galaxy spectrum templates of iKinnev et al . (1996). Since 
there are several templates provided, each having differing 
amounts of internal extinction, we again selected the tem- 
plate that produced the lowest scatter in the predicted rest- 
frame B magnitudes for L2, but between the F555W and 
F814W filters only. This optimum template had 0.15 < E(B- 
V) < 0.25 and the resulting scatter was ~0.01 mag. Again, 
no luminosity evolution correction was applied. 

Applying Equation to the HST photometry we find a 
rest-frame magnitude Mb = —20.57 + 51og/i ± 0.1 for LI, 
and Mb = —13.95 + 51og/i ± 0.1 for L2, where we use 
Ho = 100ft kms _1 Mpc _1 . We compare this to the local 
characterist ic absolute magnitu de in B measured in the rest 
frame, M t ; iRusin et alJ <l200St) use lens galaxies to find a 
value of M, — —19.7 + 5 log h ± 0.29, which implies a lumi- 
nosity for LI of 2.2 ± 0.2L», and 5.0 ± 0.5 x 10~ 3 L» for L2, 
where L* is the luminosity corresponding to M*. 

The fundamental plane is defined by 

logi? e = alogCTc +/3 < H > e +7, (2) 

in which R e is the effective radius (the radius containing 
half the total light) of the galaxy in kpc, a c is the central 
velocity dispersion of the luminous matter, < n > e is the 
galaxy's effective surface brightness in mag arcsec -2 , and a, 
(5 and 7 are empirically determined constants. The effective 
surface brightness is the mean surface brightness within the 
effective radius, and is given by 

< n > e = m (3 -^' S Q(z)+2.51og27r+51ogre-101og(l+2;), (3) 

in which r e is the effective radius in arcseconds, and tuq 
is the total magnitude of the object (i.e. the magnitude 
within an infinite aperture). We have used the total mag- 
nitude derived from the Sersic fits, on the assumption that 
the extrapolation introduces minimal additional error com- 
pared to the uncertainties in K-corrections and magnitude 
system zero points. We do not apply an evolution correction 
to < ii > e . We ado pt the values for a, f3 and 7 found by 
iBender et ail |l998) for local ellipticals, a = 1.25, /3 = 0.32, 
7 = -8.895 for H = 50 kms" 1 Mpc" 1 . The fundamental 
pla ne remains valid at redshifts comparable to that of LI 
(ee lTreu et alj ll999). Assuming that LI lies precisely on the 
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Figure 11. Velocity dispersion of LI (o"i so ) as a function of source 
redshift (z s ), assuming an isothermal mass distribution in the 
lens. The shaded area represents the velocity dispersion derived 
from the fundamental plane at the 95% confidence level when h = 
0.7. Assuming an isothermal mass distribution, we find z s > 1.2 
with 95% confidence. 



fundamental plane, we use the results from the Sersic fits to 



estimate a c 



117+?* h 



We make an assumption of isothermality for the lens galaxy 
LI in order to estimate the source redshift. Lens galaxies are 
typically found to have mass profiles within 15 % of isother- 
mal within the radius of the lensed images <| Rusin et al 



2002j : iKoopmans fc Treull200&l iRusin. Kochanek fc Keetqn 
ma 



20051: iTreu fc Koopmand 12004 iRusin fc Kochanekl 1200 A 

For an isothermal galaxy we can relate (dynamical) ve- 
locity dispersion to a function of the source redshift 
llKormann. Schneider fc Bartelmannlll994t) . because the ve- 
locity dispersion determines the angular separation of lensed 
images, A9, through 



V c J D a 



(4) 



where <Ji SO is the velocity dispersion of the isothermal mass 
(not necessarily equal to a c ), and D\ B and D B are the an- 
gular size distances from the lens to the source and from 
the observer to the source respectively. In Figure 1111 we 
have used this relation with the observed image splitting 
to plot the expected ai 30 as a function of source redshift for 
h = 0.7, assumin g that <r c and mso are approximately equal 
llKochaneklll994r>. We also plot the two source redshifts sug- 
gested bv lMcKean et al.l i2004) . which between them conve- 
niently span much of the plausible source redshift space. The 
estimated velocity dispersion of LI implies that the source 
redshift z s is greater than 1.2 at 95% confidence; a source 
with z 3 = 0.8 would require the lens to have a velocity dis- 
persion of ~330 kms -1 , assuming that LI is isothermal. 

The data available favour a high source redshift, but fur- 
ther spectroscopy will be essential if the radio source proves 
to be time-variable, since then B0631+519 would be a po- 
tential target for attempts to measure the Hubble constant. 
A spectroscopic source redshift would also be valuable for 
normalising mass models of this system. 



4 LENS MODELLING OF B0631+519 

In this Section we seek to demonstrate that the lens hypothe- 
sis is a plausible explanation for the radio emission seen from 
130631+519, by comparing the predictions of simple para- 
metric mass models to the observations. We fit mass models 
to the 5 GHz MERLIN data on Al, A2, Bl and B2 because 
these images are relatively compact. We then take the op- 
timum mass model fitted to the radio data and use it to 
invert the Einstein ring seen in the NICMO S F160W image 
by applying semi-linear inversion techniques iWarren fc Dvel 



l200a breu fc Koopman J2004 iKoopman J200rjr Finally, we 
use the mass model to examine Li's mass to light ratio. The 
unknown source redshift means that the estimated time de- 
lays, as well as Li's mass (within the Einstein radius) and 
its mass-to-light ratio are scaled by factors that depend on 
z s through combinations of the angular diameter distances 
between source, lens and observer. We therefore quote the 
affected quantities for z s = 0.80 and z s = 3.14. 



4.1 Mass modelling using the 5 GHz MERLIN 
data 

When fitting a mass model to the radio observations, the 
observed image positions are back-projected to the source 
plane through the current mass model. These back-projected 
images are used to update positions and flux densities of the 
model sources, and the updated sources are then projected 
forward onto the image plane and compared to the observed 
image positions an d flux densities t hrough a standard image- 
plane x 2 statistic <Kochaneldli99l . The statistic is defined 
over N images having angular positions 8i (i — 1...JV) and 
flux densities Si by 



N 

;=i 



e t - o' 



+ 



\si-s'i\ 2 



1=1 



(5) 



where the primed quantities are values predicted by the lens 
model being evaluated and the a P i and cr/j quantities ac- 
count for errors in positions and flux densities respectively. 
The lens model parameters are iteratively optimised by 
a do wnhill simplex method iPress. Flannerv fc Teukolskvl 
until a minimum is found. In this approach we assume 
that the observed images are unresolved, which is clearly the 
case for Al and Bl in the MERLIN 5 GHz data. A2 and B2 
are slightly resolved but are still well-represented by Gaus- 
sian models. Components Ala, Alb, Ale and Aid will be- 
come useful if future high-resolution VLBI observations are 
able to resolve Bl into corresponding images. Without ob- 
served counterparts in B, the VLBA sub-components in Al 
cannot yet be used to constrain the lens model. Images X, 
Y and Z are very extend ed and would be best treated using 
tools such as LensM EM ( Wallingtcm^Jjo^harie^fc^iaraxan 
Il996[) or LensClean dKochanek fc Naravanll992] : IWuckni~ 

boo J . 

Because Al and Bl do not arise from the same source 
as A2 and B2, this system has fewer degrees of freedom 
immediately available than for a single quadruply imaged 
source. We ignore any contribution L2 might make to the 
lens e ffect and treat LI as a singular isotherm al ellipsoid 
(SIE: iKormann. Schneider fc Bartelmannl |l994) . We allow 
the lens position to vary freely. The constraints available for 
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Model 


SIE 


SIE+shear 


Lens offset in R.A. [mas] 


+398±23 


+370 (fixed) 


Lens offset in Dec. [mas] 


— 519±5 


—520 (fixed) 


Critical radius [mas] 


603±2 


605±3 


Ellipticity 


0.13±0.03 


0.07±0.04 


Lens P.A. 


-57.1±0.1° 


-28±12° 


Shear strength 




0.03±0.01 


Shear P.A. 




-81+0.5° 


x 2 /^ 


0.41/1 


0.07/1 



Table 8. Lens models for B0631+519. Lens position offsets are 
given with respect to image Al. Position angles are measured east 
of north. The errors in the lens parameters were estimated using 
Monte Carlo simulations. 



this simple lens model are then four sets of image positions 
and flux densities. The free parameters are the positions and 
flux densities for two sources, the lens position, the lens mass 
scale and the position-angle and ellipticity of the mass dis- 
tribution. The number of degrees of freedom, v, is therefore 
1. We allow 20% errors on image flux densities to account 
for possible source variability, and we use the positions and 
position errors derived from the MERLIN 5 GHz data for 
Al, A2, Bl and B2 (because the VLBA observations resolve 
A2 and B2). The model parameters are shown in Table |H1 
and the lens is illustrated in Figure [T2l The optimum model 
has a reduced x 2 of 0.41, with roughly equal contributions 
from flux densities and image positions. 

Next, we fixed the galaxy position to that determined opti- 
cally in Section|3| and added an external shear to the model. 
This SIE+shear model was then fitted to the MERLIN im- 
age positions and flux densities. The optimum model had 
a reduced \ 2 °f only 0.07 for 1 degree of freedom, suggest- 
ing that the random errors on the constraints were over- 
estimated. Reducing the flux density errors from 20% to 
5% results in a reduced \ 2 of 0.5 which is dominated by the 
contributions from the image positions. The optimum model 
parameters do not change significantly when the flux den- 
sity errors are reduced. The optimum model parameters are 
listed in Table |H| and the model is illustrated in Figure [T2"l 
We note that the position of the best SIE-only lens is con- 
sistent with the optical lens position within the errors. The 
time delay between Al and Bl is predicted to be 10.2 h -1 
days with the SIE model or 8.5 h~ x days with the SIE+shear 
model for a source redshift of 3.14. 

The lens modelling we have performed shows that the Al, 
A2, Bl and B2 components can be understood as lensed im- 
ages of two distinct sources. It seems that Al and Bl are 
images of a compact radio core having a flat (synchrotron 
self-absorbed) spectrum, while A2 and B2 arise from a more 
extended region of emission to the west of the core, prob- 
ably a radio lobe. Images Al and A2 have positive parity, 
while images Bl and B2 have negative parity. Incorporating 
images X, Y and Z is more difficult because they are all rel- 
atively extended compared to the A and B images, and the 
radio spectra indicate that these images arise from a third 
source. Back-projecting the positions of the Y and Z images 
(from the 1.7 GHz MERLIN data) to the source plane, along 
with the brightest point on arc X, produces an approximate 
source position of +0.4l"east and — 0.49"north relative to 
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Figure 12. Lens models applied to B0631+519. The upper panel 
shows the best single SIE model. The lower panel shows the best 
SIE+shear model when the lens galaxy position is fixed to the 
measured optical/infrared galaxy position. The ellipse in each 
panel is the 4-image critical curve. The dashed contour lines rep- 
resent the time delay surface of the lens relative to image Al. The 
contours in the upper panel begin at 1.23 h^ 1 days (for the first 
contour surrounding A1/A2) and increase in steps of 2.05 ft. -1 
days. In the lower panel, the first contour surrounding A1/A2 rep- 
resents 1.11 hr 1 days and the contours increase in steps of 1.88 
h~ 1 days. The contour levels are calculated for a source redshift 
of 3.14, and should be multiplied by a factor of 3.5 for z s = 0.80. 
The lensed images are represented by filled circles. The inset di- 
agram shows the source plane; the astroid curve is the 4-image 
caustic, and filled circles represent the modelled source positions. 
The estimated position of the source responsible for images Y, Z 
and the arc X is marked by a light grey filled circle. 



Al (marked on Figure ^J. The model thus suggests that 
the X, Y and Z images are produced by a counter-lobe on 
the eastern side of the source that happens to lie inside the 
4-image caustic. With only the VLA data, showing X and 
Y but no Z, a model involving emission only on the western 
side of the source's core might also be plausible. 
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4.2 Non-parametric de-lensing of the Einstein 
ring 

The blue color of galaxy L2 and its consequent apparent 
absence in the HST-NICMOS F160W image, allows us to 
make use of the infrared Einstein ring to further test the 
validity of the lens mass model of CLASS 0631+519, which 
is based on the MERLIN radio data alone. 

The Einstein ring is de-lensed (see IWarren fc Dvd l2003t 
iTreu k. Koopmansl 12004 iKoopmansl l200ElL using"the SIE 
plus shear model from Table |H] It is stressed that the lens 
mass model is not adjusted in any way to better fit the ob- 
served Einstein ring. A Tiny Tim generated PSF is used to 
account for blurring of the image. The structure of the Ein- 
stein ring is remarkably well reproduced (Figure USTl . given 
that the lens mass model was only fitted to the radio jet 
structure. Despite the overall good quality of the recon- 
struction, a clear deviation from the model is found in the 
residual image, just north of LI on the Einstein ring and 
slightly offset from the position of L2. Neither varying the 
mass model, nor alleviating the regularization can remove 
the feature completely. 

We propose two possible explanations for this apparent 
"anomaly" : (a) it is due to residual emission from L2 - still 
present at a low level in the NICMOS-F160W image - or a 
dwarf-galaxy associated with LI, or (b) it is due to a small- 
scale perturbation of the lens mass m odel by the prese nce of 
L2, creating an anomaly in the ring jKoopmansl l 2005 ) sim- 
ilar to flux-ratio anoma lies in compact lensed sources (e.g. 
iMao fc Schneideilll998f) . Although the first explanation (s) 
can not be completely excluded, the anomaly is situated 
near the brightness peak of the Einstein ring, appears to fol- 
low the curvature of the ring and is offset from the centroid 
of L2. Each of these observat ions seem to supp ort the sec- 
ond explanation. As shown in IKoopmansl i2005t) . one would 
indeed expect such anomalies in Einstein rings, as a result 
of tiny perturbations of the dominant lens potential. The 
foreground galaxy L2 could provide the required perturba- 
tion. A full analysis of the nature of the anomaly, however, 
is beyond the scope of this paper. 

4.3 Mass to light ratio of LI 

We have sufficient information to estimate the B-band mass- 
to-light ratio of LI in the region enclosed by the lensed im- 
ages. The mass M between the lensed images of a source 
(in the case of a circularly symmetric lens) is related to the 
image separation A9 and the angular diameter distances Di 
(observer to lens), D s (observer to source) and Di s (lens to 
source) by 

DiD a 



CLASS B0631+519 



M = 3.074 x 10 7 



(Aey 



(6) 



where A6 is given in arcseconds. The estimate is signifi- 
cantly affected by the unknown source redshift (but not 
by assumptions concerning the distribution of mass inside 
the radius of the images). The ratio of angular diame- 
ter distances in Equation [(J decreases by a factor of ~3.5 
when the source redshift is increased from 0.8 to 3.14. 
We therefore estimate the lensing mass usi ng both of the 
sugges ted redshifts for the source listed in iMcKean et alJ 
(2004), since they conveniently span most of the range of 
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Figure 13. Non-parametric lens reconstruction of CLASS 
B0631+519. Top-left: The galaxy-subtracted NICMOS image of 
the Einstein ring. Top-right: The best reconstruction of the ob- 
served Einstein ring. Bottom-left: The residuals between obser- 
vations and the model. Bottom-right: The reconstructed source 
brightness distribution. For the source reconstruction, we use a 
30 X 30 pixel grid with a scale of 0.02 arcsec pix — 1 . We mask 
all image pixels below 3 cr (indicated in the image reconstruction; 
a = 0.011 counts pix - 1 ). The regularization parameter for the 
source is A s = 0.04, resulting in x 2 /d.o.f. ~ 1. We use a regu- 
lariza tion matrix that minimizes the curvature of the source (see 
lKoopmansll2005l) . Note that a significant "anomaly" is left in the 
residual image (see text). The two circles in the residual image 
indicate the centroids of LI (at the origin of the plot) and L2 at 
—0.33 arcsec east, +0.74 arcsec north of LI 



plausible source redshifts. The separation between images 
Al and Bl is 1.162+0.003 arcsec, which at the redshift of 
LI corresponds to a linear distance of 5.6/i _1 kpc. From 
EquationHJwe estimate the mass of LI within the images at 
2.17+0.01 xlO 11 /! -1 M @ using a source redshift ofz a = 0.80 
or 6.27 ± 0.03 x 10 10 /i _1 M Q using z 3 = 3.14. Combining 
these estimates with the rest frame B-band absolute mag- 
nitudes derived from the Sersic profile fits in Section 13.31 
(corrected for a 1.16 arcsec diameter aperture) produces es- 
timates for the mass-to-light ratio of LI in the B filter of 
15±2/i Mq/Lq using z 3 = 0.80, or 4.3 ± 0.4ft Mq/Lq using 
z 3 = 3.14. We have taken the absolute magnitude of the Sun 
in B to be 5.54. The mass-to-light ratio estimates apply to 
the central region of LI encircled by the lensed images (i.e. 
that within a mean radius of 2.8 ft -1 kpc from Li's centre). 
Both figures are comparable with M/L estimates from lo- 
cal ellipticals jGerhard et al.ll200lf) . although the figure for 
z a = 3.14 is relatively low given h = 0.7. 



5 CONCLUSIONS 

Together the radio and optical/near-infrared observations 
presented here show that CLASS B0631+519 can be under- 
stood as a gravitational lens system. A simple SIE+shear 
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mass model successfully reproduces features of the radio 
observations, and a non-parametric approach to the near- 
infrared data has identified the perturbing effect of L2 on 
the smooth lensing properties of the system. The mass mod- 
elling done so far leads us to conclude that the background 
radio source probably consists of three radio components 
in a classic double-lobed core-jet arrangement. The eastern 
lobe is located inside the quad-image caustic for the system 
and is lensed into the arc X and the more compact images 
Y and Z. The core region is doubly imaged into Al and Bl, 
while the western lobe is lensed into images A2 and B2. The 
magnification provided by the lens and the high resolution 
of the VLBA has allowed identification of the flat-spectrum 
AGN core and revealed jet components from both the ap- 
proaching and the receding jets. The near-infrared source 
centre coincides with the radio source positions to within 
the errors. 

The wealth of constraints available means that this sys- 
tem will benefit from application of deconvolution/inversion 
techniques that can utilise the extended images seen in the 
radio data. Having a good mass model for this system will 
allow us insights into the mass distributions of ellipticals at 
redshifts of 0.6, information that is not readily available in 
any other way. The hints of radio source variability mean 
that it may be possible to determine a radio time delay and 
use it to measure the Hubble constant. 

This lens system is relatively unusual in having two galaxies 
of different redshifts present along the line of sight to the 
source, a situation that also occurs in CLASS B2114+022 
llKing et alJll99flt lAueusto et alJ l20Qlft an d in PKS 1830- 
211 iPramesh Rao fc Subrahmanvanlll988h . In B21 14+022 
the effect of th e nearer lens is much more significant than 
in B0631+519 llChae. Mao fc Augustdliool . In PKS 1830- 
211 the nearer lens was detected through its HI absorption 
jLovell et al.l |l996t) at a redshift of 0.19 (the main lens- 
ing galaxy has a redshift of 0.886). The correct interpre- 
tation of the optical data in PKS 1830-211 is presently un- 
clear, although the main lens is certa inly a face-on spiral 
JWinn et alJl2002l; ICourbin et al"ll2002ft . 

Although L2 appears to have a relatively minor influence on 
the lensing properties of B0631+519, the radio emission from 
the lensed source must pass through L2's ISM. Searches for 
molecular absorption lines may be worthwhile to provide 
extra information on the nature of this galaxy. The mass 
distribution in L2 may well be constrainable with further 
non-parametric modelling work. 
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